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Abstract. The SIM Lite mission will undertake several planet surveys. One 
of them, the Deep Planet Survey, is designed to detect Earth-mass exoplanets 
in the habitable zones of nearby main sequence stars. A double blind study 
has been conducted to assess the capability of SIM to detect such small planets 
in a multi-planet system where several giant planets might be present. One 
of the tools which helped in deciding if the detected planets were actual was 
an orbit integrator using the publicly available HNBody code so that the orbit 
solutions could be analyzed in terms of temporal stability over many orbits. In 
this contribution, we describe the implementation of this integrator and analyze 
the different blind test solutions. We discuss also the usefulness of this method 
given that some planets might be not detected but still affect the overall stability 
of the system. 



A d ouble blind study h as been conducted in the framework of the SIM-Lite 
mission QTraub et al||2009[ see also Traub et al. in this volume) to assess the 
capability of SIM to detect such small planets in a multi-planet system where 
several giant planets might be present. Five independent teams of dynamics ex- 
perts generated ensembles of hundreds of plausible planetary systems that might 
form around solar-type stars. A data simulation team generated realistic simu- 
lated astrometric and radial velocity measurement data sets for 60 stars chosen 
from the SIM target list, by perturbing each with a planetary system randomly 
drawn from the ensembles generated by the dynamics teams, and then adding in- 
strument noise. Next, four independent analysis teams processed the simulated 
data sets, detecting and fitting the orbits and masses of the planets orbiting each 
of the 60 stars. Finally, the analysis teams submitted their consensus solution 
for each star, chosen after cross-comparison and discussion of the solutions, and 
before being allowed to see the true solutions. 

We used HNBody, a symplectic in tegration package for hierarchical N-Body 
systems (version 1.0.3) developed by |Rauch fc Hamilton| ( [2002] ) . It integrates 
the motion of particles in self-gravitating systems where the total mass is dom- 
inated by a single object; it is based on symplectic integration techniques in 
which two-body Keplerian motion is integrated exactly. HNBody is primarily 
designed for systems with one massive central object and has been u sed previ- 
ously for extrasolar planet simulations ( |Veras fc Armita ge 2005, 2006[ ). We used 
the parameters given by the different teams, namely the mass of the planets, 
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Figure 1. Planetary system #13 of batch 2/phase 2: xy-, xz- and yz- 
projections of the orbits for the solutions of the 4 teams. 

their period translated in semi-major axis using the Keplerian formula, the orbit 
eccentricity, the orbit inclination, the longitude of the ascending node, the argu- 
ment of periapse in years, and for the periapse passage time, since HNbody only 
accepts the epoch and time, we rounded the time of periapse to get the Epoch 
and the time was then just the remaining decimals. In order to ensure that we 
sample correctly the orbit evolution, the total integration time has been set to 1 
million years while the elementary step has been set up to 1/20 of the shortest 
period of the system. HNbody computes then the relative errors in energy and 
angular momentum exhibited by the integration. If the error mean and rms are 
less than 10 -4 over the total integration time, the system was considered stable. 

The plots presented in this Fig. [T] is representative of the temporal evolution 
of the orbit parameters found by the different teams and their 3D orbits. Figure 
[T] display the system for which there was intense discussion about the different 
solutions and for which the stability would be decisive in choosing the solution. 
This system was found stable (energy and momentum errors less than 10 -8 
in 10 6 yrs) even though there were secular changes eccentricity for the closest 
planet. 

The question arises whether such a stability test is appropriate for validating 
a solution. Obviously, if a system misses large and massive planets orbiting 
further away, then a stable system might be found unstable because incomplete. 
This is the case of the system #3 in batch #3 of phase 2. The unraveled solution 
is in fact a 9-planet system which is stable when HNBody integrates the orbit over 
1 million years. The reverse is true especially if only one planet is discovered, 
because integration of only one planet is very likely to be stable. Nevertheless, 
the orbit integration tool appeared to be useful to choose the final solutions. 
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